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METHOD FOR DETERMINING THE HUMIDITY AND DENSITY OF A 
DIELECTRIC MATERIAL 

Background of the Invention 

The invention relates to a method for determining the humidity and/or density 
of a dielectric material in a resonator filled with the material, comprising a 
sender and a receiver. 

The dielectric properties of a material, described by the complex relative 
dielectric constant £r = s[- iCf ^ can be affected by humidity and density in case 
of porous materials. Humidity and density thus modify the scalar parameters 
resonant fi-equency fr and resonator quality Q of a material-filled resonator in 

contrast to those of the empty, air-filled resonator (fro Qq) in such a way 
that 

f _ ^^0 1 1 fir 

~ and Q^IT = q7 IT, wherein 1 < e[ < ej^ and 0 ^ e" ^ CrL, wherein 

and max are the maximum values resulting for the respective material 
from the assigned humidity range and density range. 

In prior art, different methods are known for measuring the humidity or density 
of granular materials, in which the resonance behavior of a matter-filled 
resonator is used. 

For example, from US Patent No. 5,666,061 a method is known for measuring 
the humidity in granular materials by means of microwaves. There, one 
electronically reacts to a threshold while bidirectionally sweeping a frequency 
range. From the chronological properties of pulses created in such a way the 
components of the dielectric constant are deduced. 
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All methods are relatively slow and have either a dependency of the density of 
the material or have relatively large errors in the determination of the humidity. 

The invention is underlied by the problem to specify a method of the type 
initially mentioned, by which a fast, accurate and density-independent 
determination of the humidity is possible. 

Sunraiary of the Invention 

For the purpose of the invention, any arbitrary quantity describing the 
resonance width can be considered as a bandwith in the following, wherein 
appropriate adjustments have to be provided in the respective definitions, in 
particular of the threshold values, and equations. 

By digitally recording the resonance curve, a fast acquisition of the measuring 
values is possible. Thus, a close chronological-spatial assignment of the 
humidity determined from these values to material dynamically guided through 
the resonator is possible. For this purpose, the invention provides for that the 
sender emits a signal; a resonance curve of the filled resonator is sweeped, 
wherein respective related signal strength values of the receiver signal are 
measured at different fi-equencies; the resonant fi-equency and the bandwidth 
are determined for the filled resonator fi*om the measured points; and the 
humidity and/or density of the material is calculated by solving a second 
system of equations comprising the resonant fi-equencies and the bandwidths of 
the empty and of the filled resonator and known calibration coefficients of the 
resonator. This method enables a density-independent determination of the 
humidity of the material. Besides, the density of the material can be acquired 
with httle effort. 
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A preferred embodiment provides for that, from the points for determining the 
bandwith of the filled resonator, either the quantities resonant frequency, 
resonator quality and resonance maximum are determined and the bandwidth is 
calculated therefrom, or cut-off frequencies are determined and the resonant 
frequency and the bandwidth are calculated therefrom. 

In another embodiment a lower threshold value is calculated and a second 
sweeping pass with smaller step sizes is performed in that range in which the 
signal strength values are higher than the threshold value. By a two-pass 
procedure the accuracy can be significantly increased and, nevertheless, the 
required time can be kept short, in particular if, during the second sweeping 
pass, a decreased step size is used in a range of the resonance peak only. 

Advantageously, sweeping the resonance curve is performed in equally spaced 
steps. The simplest and fastest form of sweeping consists of equally spaced 
steps. Variable step sizes in the second sweeping pass can shorten its duration. 

Preferably, the sender is operated using a constant strength. The measuring 
values at the receiver can be used without adjusting or scaling them in case of a 
constant signal strength. 

In a possible embodiment, the cut-off frequencies of the resonator are 
determined by determining the point having the highest receiver signal strength 
value, and, starting from this point, calculating a threshold value; determining 
two respective proximate points for positive and negative slope sections, the 
signal strength values of these points lying below and above the threshold 
value; calculating first and second cut-off frequencies therefrom by 
respectively interpolating between the proximate points. Determining the cut- 
off frequencies by interpolating between point pairs surrounding a threshold 
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value as an initial parameter for determining the humidity and/or density is a 
fast and simple method. 

An advantageous embodiment thereby provides for that the threshold value 
corresponds to an attenuation of 3 dB in relation to the highest signal value. If 
the threshold value is chosen corresponding to an attenuation of 3 dB, starting 
from the maximum signal strength value, the equations to be solved obtain a 
very simple form. 

In an alternative embodiment the quantities resonant frequency, resonator 
quality and resonance maximum of the resonator are determined by arbitrarily 
and/or randomly selecting three points and solving a first system of equations 
for these quantities, the system consisting of three equations of an analytic 
resonance curve that are valid for the three points. The resonant frequency, the 
resonator quality and the resonance ampitude can be directly determined as the 
initial parameters for determining the humidity and/or density even faster and 
with less error than by interpolation by solving a first system of equations 
consisting of three equations and thus being completely determined, wherein 
for each point out of a group of three one of the three equations is valid, the 
three points being selected from the present set. In particular, noisy resonance 
curves can be analyzed more exactly this way. 

In another alternative embodiment the quantities resonant frequency, resonator 
quality, resonance maximum of the resonator are determined by arbitrarily 
and/or randomly selecting a set of points whose number is an integer multiple 

of three and at least six, and splitting up the point set into three equally sized 
groups; for each combination of three points, wherein each point comes from a 
different group, solving a first system of equations for these quantities, the 
system consisting of three equations of the analytic resonance curve valid for 
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these three points; and creating the average for each quantity from the values 
calculated at the combinations. Even more exact values than with three points 
can be obtained by creating several triple-groups from the present point set and 
averaging the initial parameters obtained thereby over all groups. 

In an advantageous embodiment, as a condition for arbitrarily and/or randomly 
selecting the points, the signal value of a point to be selected is higher than the 
highest signal value attenuated by 3 dB. In both procedures using three points 
or a multiple thereof, preferably only those points from the resonance curve are 
selected whose signal strength is higher than the maximum value of all 
measurement values, attenuated by 3 dB, as this way only significant values are 
used. 

The system of equations to be solved is preferably chosen in such a way that it 
describes, in a good approximation, the correlation of humidity and density 
with a variation of resonant frequency and resonator quality or, respectively, 
with a variation of resonant frequency and bandwidth. 

Advantageously the second system of equations is non-linear. 

In a preferred embodiment the sweeping by means of the sender is performed 
up to the microwave area. 

Advantageously, voltage values or current values of the receiver are used for 
measuring the receiver signal. Preferably, the electrical voltage rising at the 
receiver serves for acquiring the resonance curve, because it is measurable 
easily and without back coupling. However, the current in a receiver circuit can 
be measured, too. 
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In the following, the denomination Ui is used for the measured quantities. 
However, in doing so, it shall not be implied to use the voltage only. 

The resonance curve to be reconstructed from the discrete measuring points is 
an approximation as-good-as-possible to the real resonance curve only. It is 

characterized by the course of the signal amplitude U at the receiver against the 

supplied frequency ^, wherein Ur is the resonance maximum and Q is the 
quality of the resonance: 

u= , 

The cut-off frequencies fa and are those frequencies at which the signal 
strength exceeds and falls below a defined value, respectively. For this purpose, 
preferably the value corresponding to an attenuation of the maximum by 3 dB 
U = U — 

is chosen: ^ " " /2 . Between the cut-off frequencies and the curve 
parameters resonant frequency fr and resonator quality Q there are the 

correlations: ^ Vfa^ , Q = fj-fa . The distance between the cut-off 
frequencies is defined as the bandwidth = ft - fa = q* 



In the following, the invention is described in further detail using examples of 
embodiments. 

Brief Description of the Drawings 

Figure 1 shows a sweeped resonance curve; 

Figure 2 shows a resonance curve sweeped in a two-pass procedure; 
Figure 3 shows a first way for determining the resonance parameters; 
Figure 4 shows a second way for determining the resonance parameters; 
Figure 5 shows a third way for determining the resonance parameters; 
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and 

Figure 6 shows a two-dimensional representation of calibration curves. 
Detailed Description of the Invention 

In Figure 1, two resonance curves are depicted above each other like they have 
been recorded using the method according to the invention; the right resonance 
curve at an empty, the left resonance curve at a material-filled resonator. 

By sweeping (wobbling) the resonator between the start frequency fstaiti and 
the stop frequency fstopi, the resonance curve is acquired in discrete steps. The 
start frequency is calculated from the maximally shifted frequency 

*^ V^^max ^ the maximally changed quality OlIT ^ ^r^ax and the normalized 

voltage ratio ^'^^^ Jl : ^stanl - 2Q„ V ^ iyiQj ' +Ir„ jYiQ 

stop frequency, using fr = fro, results as: 

f _ / , // fio \7 . eo' 

Islopl --2Q„, V a2 '^V'^2QniJ • a2 + ^ro 

The sweeping speed depends, among others, on the number of sweeping 
points and can be increased by a two-pass procedure as depicted in Figure 2. 
For this purpose, first a sweeping is performed using a smaller number of 

sweeping points and, hence, a larger frequency step size - n,-i and, in 
a second iteration, after previously determining the start and stop frequencies 

fstara =/| U^Va A f>Q, fstopl = f\ V^Ua A/<0 or fstara =/| V^Ua A U<Ur^ 

fstopi-f\u^UaMJ>Ur^ another sweeping is subsequently performed using a 

A r» fstop2~fstart2 

smaller frequency step size - nj-i between these frequencies , 
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From the sweeped, measured resonance curve the resonator parameters 

can be determined, for example, according to the fafb procedure, 
the three-points procedxire or the 3k-points procedure. Other procedures are 
possible, too. 

Figure 3 shows the fafb procedure. It is based on directly determining a first 
and a second cut-off frequency fa> fb from the measured resonance curve. For 
this purpose, first the point (also called sweeping point) 
(f = fiiiax/U = Uniax)having the highest voltage U^Umax is determined for 

IT _ ^max 

calculating a 3 dB threshold line 72 . Subsequently, the respective two 

points in the immediate proximity of the 3 dB threshold line are determined. By 
linearly interpolating between the sweeping points al and a2 as well as bl and 
hi the first and the second 3 dB cut-off frequency fa, fb are obtained: 

fa = fal + Ua2-Ua, (fa2 " fal)^ fb = fbl + Ub2-Ubi (^^2 " fbl). 

From the equations ~ V^Wb , Q = fj-fa ^ fr and Q can be calculated. The 
resonance voltage then results as: 

Ur = U™«-^l+Q<^-i)' 

Calculating the cut-off frequencies is more or less error-prone because of the 
discretization of the resonance curve - on one hand, by the determination of the 
3 dB threshold line (Ua = Umax/ 4l ^ Ur/ 42 ) derived from Umax and, on the 
other hand, by the linear mterpolation between the sweeping points al, a2, bl, 
b2. A higher sweeping rate indeed reduces the errors, but also increases the 
sweeping time and thus reduces the sweeping speed. At best, a non-linear 
interpolation is possible in this procedure in order to improve accuracy. 
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In the three-points procedure shown in Figure 4, three sweeping points 

U,,U2,U3 

are selected arbitrarily or by a random generator preferably above 
— 

the line Us = S' Umax with 7?. 



By solving the first system of equations Gl: 

the resonator parameters/*-* Q ^"'^ ^roie obtained: 

_ l u?(u^ui)ff+ul(u?-u^)fl+u|(ui-u?)ff 

V fi ft 



Q = 



Ur=U,^l+Q<f-f)' 



The advantage of the three-points procedure in relation to the fafb procedure 
lies in that no errors occur because of the discretization of the resonance curve, 
and that only few points are needed for sweeping the resonance curve, and thus 
a high sweeping speed is achieved. 

If the resonance curve is noisy, error-prone resonator parameters are 
determined by the three-points procedure, too. In order to reduce the impact of 
the noise, more than three sweeping points are selected, namely a 
number^-k, k = 2,3,4,...^ and are split up into three groups Ml, M2, M3, 
having k points per group above the line Us = s • Umax . Figure 5 shows this 3k- 
points procedure. For all k^ possibilities of combining one respective point 
from each of the three groups, the resonator parameters are determined from 
the first system of equations Gl according to the three-points procedure 
described above, and the k^ partial results are averaged subsequently: 
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fr = F 2?=, f, . Q = 1^ l£ Qj . U, = 1^ U, . 

In order to determine humidity and density from determined resonator 
parameters, the used arrangement is caUbrated using material of known 
humidity and density. The calibration is advantageously performed out of the 
regular operation. 

hi the calibration, the electrical resonator values resonant frequency fr, 

resonator quality Q and bandwidth BW(BW = fr/Q) are respectively assigned 

to the material quantities humidity ¥ and density Q. This assignment between 
the material values and the electrical values is non-linear and can be described 
for a defmed range of humidity and density with sufficient accuracy, for 
example, by the following second system of equations G2: 

Afr = HiQ^W^ + ^UQW^'^HiQ^W'^HxQV^^faQ^-^^frxQ 

ABW = aewz^'^ + aew, Qy/'^^ h^^^g^ y/ + bsw, Qy/^ cbwzQ^ + cbw, Q 

Afr and ABW are the differences of resonant frequency and bandwidth, 
respectively, between the empty arid the material-filled resonator: 

Afr = fro - fr., ABW = BWn, - BWq 

Now, it is an object of the calibration to determine the twelve calibration 
coefficients ^fn^ afr2» bf,,, bf,2» Cfru ^fa^ ^bw,, aawj, bew,, bewj, csw,, cewj 
for the respective material from a sufficient large number of calibration values 
(Afr and ABW ^j^h related material values <^ und 6), 

For this purpose, the measuring values of the resonant frequency difference Afr 
and the bandwidth difference ABW are assigned to the humidity and density 
values that have been determined by appropriate reference methods. These 
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calibration values characterizing the resonant frequency and bandwidth 
depending on humidity and density are the basis for determining the calibration 
coefficients. The regression of the calibration values of equal humidity is 
performed in the form of Afr = df,2e^ + df,ig^ and Afr = dBW26>^ + dBW,6>, 
wherein the regression curves have to run through the origin, as and ABW 
are zero for the empty resonator (air having Q - 0). The regression provides 
calibration points of equal density depending on the humidity. Using these, 
another regression of the form Afr= af,^^^ +bf,^ + Cf, and 
ABW= aawV^^+bBw^^^ + CBw is performed. The values determined this way 
for the regression coefficients ^fr? bf,, Cf, and ^bw, baw, Cbw are plotted 
against the density and, therefrom, the calibration coefficients are determined 
by a quadratic regression. 

f4g»-Figure_6 shows the results of the calibratrion in the form of a set of 
curves which also serves for analysis. 

From the measiu-ed resonant frequency difference values Afr and bandwidth 
difference values AB W humidity content ¥ and density Q are calculated for the 
respective material by solving the above second system of equations G2. 

In doing so, two real and two imaginary roots [[do]] resuh. From the course of 
the calibration curves in the Afr-ABW diagram it can be determined if there is 
only one real solution in the humidity and density range of interest. In the Afr- 
ABW diagram the bandwidth difference is represented in dependency of the 
resonant frequency difference for curves of equal density and humidity. If the 
course of these curves, characterized by the points A, B, C und D in Figure 6, is 
continuous and unique in the humidity and density range of interest then only 
one real solution exists in this range. 
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For solving the second system of equations G2 an iterative procedure is 
appropriate. For this purpose, the second system of equations G2 is solved for 

_ _ bfr2(H-bfri W bfr2g+bfri ^ _ Cfr2g^-K:fr|-Afro 

_ ^BW2g+bBW| \f bBW2^bBWi ^^ CbW2('^+CBW,-ABWo 

^ 2aBW2(?+2aBW]^ ^2aBW2(?+2aBW,(? ) ~ aBw^^^^i^Bwji 

From the intersection of both equations in a humidity-density diagram die 
sought values for and Q are obtained. 

10 
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List of reference numerals 

Resonant frequency in the general case 

fro Resonant frequency of the empty resonator 

5 fr™ Resonant frequency of the filled resonator 

fstarti , fstopi Start and stop frequencies for first sweeping pass 

^max Highest signal strength value measured 

a, s Threshold value factors 

fstart2, fstop2 Start and stop frequencies for second sweeping pass 

10 fmax Frequency at which the highest signal strength 

value is present 

(fai/Uai), (fa2/Ua2) Proximatc points of the first cut-off frequency 

(fbi/Ubi), (fb2/Ub2) Proximate points of the second cut-off frequency 

(f i/U i) ... (fs/Ua) Three selected points 

15 Ml, M2, M3 point groups 

(fii/Uu) ... (fi4/Ui4) Elements of point group Ml 

(f2i/U2i) ... (f24/U24) Elements of point group M2 

(f3i/U3i) ... (f34/U34) Elements of point group M3 

A,B,C,D Limits of the hvunidty and density range of interest 
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